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Magnetic and crystallographic properties of the mineral langite Cu 4 ( 0 H) 6 S 04 - 2 H 20 are reported. 

Its layered crystal structure features a peculiar spatial arrangement of spin-| Cu^^ ions that 
arises from a combination of corner- and edge-sharing chains. Experimentally, langite orders 
antiferromagnetically at Tjv — 5.7 K as revealed by magnetization and specific heat measurements. 

Despite this very low energy scale of the magnetic transition, langite features significantly stronger 
couplings on the order of 50-70 K. Half of the Cu^'*' spins are weakly coupled and saturate 
around 12 T, where the magnetization reaches 0.5/rs/Cu. These findings are rationalized by 
density-fnnctional band-structure calculations suggesting a complex interplay of frustrated exchange 
conplings in the magnetic planes. A simplified model of coupled magnetic sublattices explains the 
experimental features qualitatively. To start from reliable structural data, the crystal structure of 
langite in the 100-280 K temperature range has been determined by single-crystal x-ray diffraction, 
and the hydrogen positions were refined computationally. 

PACS numbers: 75.50.Ee,75.10.Jm,71.20.Ps,61.50.Ks 


I. INTRODUCTION 

Low-dimensional magnets show unique diversity of 
crystal structures and associated spin lattices, where a 
plethora of quantum phenomena can be observed. The 
physics of quantum spin chains has been actively ex¬ 
plored in Cu^“'" compounds featuring chains of corner- or 
edge-sharing CUO4 plaquette units. The corner-sharing 
geometry results in uniform spin chains with a negligi¬ 
bly small second-neighbor coupling, as in Sr2Cu03,^’® 
AgCuV04® and KCuMo04(OH).^ The edge-sharing ge¬ 
ometry is by far more common. It gives rise to com¬ 
peting nearest-neighbor and next-nearest-neighbor cou¬ 
plings, where the former (Ji) is typically ferromagnetic, 
while the latter (J2) is antiferromagnetic. Such Ji — J2 
frustrated spin chains develop incommensurate spin cor¬ 
relations and helical magnetic order,®’® although few in¬ 
stances of ferromagnetic intrachain spin order are known 
as well.^®’^^ The helical spin arrangement observed in 
simple binary compounds CuCl2^® and CuBr2^® and in 
more complex materials like linarite PbCu(0H)2S04,^^ 
all being frustrated Ji — J2 spin chains, may trigger elec¬ 
tric polarization induced by the magnetic order, thus 
leading to multiferroic behavior.Additionally, small 
interactions beyond the isotropic Heisenberg model lead 
to an intricate magnetic phase diagram, including mul¬ 
tipolar (three-magnon) phases, which has been studied 
recently.^® However, the complex interplay of frustration 
and anisotropy needs further investigations on different 
systems as, e.g., LiCuV04.^®’®® 

One may naturally ask what happens when two types 


of spin chains, those with edge- and corner-sharing ge¬ 
ometries, are placed next to each other within one 
material. Spin systems comprising several magnetic 
sublattices with different dimensionalities and energy 
scales may have very unusual low-temperature proper¬ 
ties. When two sublattices are weakly coupled, they 
are, to a certain extent, independent, hence two mag¬ 
netic transitions manifesting the ordering within each 
of the sublattices could be observed. On the other 
hand, the ordering within one sublattice will neces¬ 
sarily depend on the other sublattice, because three- 
dimensional ( 3 D) long-range order typically involves in¬ 
teractions between the sublattices. Unusual manifesta¬ 
tions of quantum order-from-disorder have been observed 
in Sr2Cu304Cl2®^^®® featuring interpenetrating square 
lattices with drastically different exchange couplings. In 
CUP2O6, where spins, arranged on a planar square lat¬ 
tice, coexist with uniform spin chains, very strong spin 
fluctuations are observed even below the Neel temper¬ 
ature T/v, and the value of Tpf is unusually low for a 
quasi-two-dimensional ( 2 D) antiferromagnet.®® The co¬ 
existence of corner- and edge-sharing 00®+ chains could 
be even more interesting because of the different nature 
of spin correlations, which are expected to be antiferro¬ 
magnetic collinear and helical, for the corner- and edge¬ 
sharing chains, respectively. 

The respective magnetic ground state of these com¬ 
pounds depends very subtly on the interplay of vari¬ 
ous exchange integrals, including possible frustration and 
strong quantum fluctuations. In particular, in edge¬ 
sharing geometries (with Cu-O-Cu bond angles near 
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90°) the leading exchange integrals and, thus, the actual 
magnetic model are often difficult to establish due to a 
pronounced dependency of the exchange on the struc¬ 
tural details: Small changes of bond angles or minor 
changes of the local Cu-0 environment, e.g. by attached 
H-atoms,^^’^® may even swap the ground state qualita¬ 
tively. Owing to the high complexity of the structure- 
properties relation in these compounds, the combination 
of experimental investigations with theoretical methods 
appeared to be very successful to disentangle the com¬ 
plicated interplay. In particular, in recent years den¬ 
sity functional calculations have developed to a valuable 
tool, establishing accurate magnetic models on micro¬ 
scopic grounds,even for involved geometries like 
coupled edge-sharing chains or various magnetic sublat¬ 
tices in a single compound. 

The coexistence of the edge- and corner-sharing 
geometries is rather common for Cu^“'" minerals. In 
antlerite^^ and szenicsite,^^ one edge-sharing chain 
is encompassed by two corner-sharing chains that 
together form a three-leg spin ladder. The edge- 
and corner-sharing chains can also form infinite 
layers, as in deloryite Cu 4 (U 02 )(Mo 04 ) 2 (OH) 6 ,^^ 
derriksite Cu 4 (U 02 )(Se 03 ) 2 ( 0 H) 6 ,^^ niedermayrite 
Cu 4 Cd(S 04 ) 2 ( 0 H) 6 ,^^ and langite Cu 4 ( 0 H) 6 S 04 - 21120 . 
The crystal structure of the latter mineral is shown in 
Fig. I. Layers of alternating edge- and corner-sharing 
chains (Fig. 1, right) are well separated by SO 4 sulphate 
groups and water molecules. A somewhat similar 
structure without water molecules has been reported 
for the mineral brochantite Cu 4 ( 0 H) 6 S 04 ^®’^^ that, 
however, features a much smaller interlayer separation, 
hence substantial interlayer couplings can be expected. 
In this paper, we focus on the magnetism of langite, 
where individual structural planes should be very weakly 
coupled magnetically and sufficiently pure natural 
samples of this mineral are available. 

We also show that, in contrast to our expectations, in¬ 
dividual structural chains in this material cannot be con¬ 
sidered as weakly coupled magnetic chains. On the other 
hand, two sublattices formed by the structural chains of 
either type, can be distinguished in the overall very in¬ 
tricate spin lattice. These sublattices reveal drastically 
different magnetic couplings and facilitate the description 
of the complex low-temperature magnetism on a qualita¬ 
tive microscopic level. 

The paper is organized as follows: Applied experimen¬ 
tal and theoretical methods are described in Sec. II. The 
crystal structure of langite including the single-crystal 
data collected at low temperatures and hydrogen posi¬ 
tions determined computationally is presented in Sec. III. 
Sec. IV provides experimental results on thermodynamic 
properties of langite. The electronic band structure and 
computed exchange coupling constants are discussed in 
Sec. V. Eventually, a detailed discussion and summary 
are given in Secs. VI and VII, respectively. 


II. METHODS 

All experiments have been performed on a natural sam¬ 
ple (Fig. I) of langite from the Podlipa and Reinera mine, 
Lubietova, Slovakia. The sample quality was first thor¬ 
oughly controlled by laboratory powder X-ray diffraction 
(XRD) (Huber G670 Guinier camera, GuKq, 1 radiation, 
ImagePlate detector, 29 = 3 — 100° angle range). 

Single-crystal X-ray diffraction between 100 and 280 K 
was performed on a Bruker SMART APEX CCD- 
diffractometer equipped with a Cryosteam liquid nitro¬ 
gen low-temperature device. A single crystal, selected on 
the basis of its optical properties (sharp extinctions, reg¬ 
ular shape and homogeneity in color) was glued on top of 
a glass capillary (0.1 mm). Intensity data were collected 
with graphite-monochromatized MoKq, radiation (50 kV, 
30 mA). The crystal-to-detector distance was 40 mm and 
the detector was positioned at —28° 20 using an w-scan 
mode strategy at four different (^-positions (0°, 90°, 180° 
and 270°). 630 frames with Aw = 0.3° were acquired for 
each run. The 3D data were integrated and corrected 
for Lorentz polarization and background effects using 
the APEX2 software (Bruker - Nonius, 2004). Struc¬ 
ture solution (using Patterson methods) and subsequent 
weighted full-matrix least-square refinements on were 
done with SHELXL-97^® as implemented in the program 
suite WinGX 1.64.39 

All further experiments were performed on a powder 
sample since the natural crystals are very small and frag¬ 
ile as well as strongly intergrown, preventing us from col¬ 
lecting a sufficient amount of single crystallites for mag¬ 
netic and specific heat measurements. Magnetization 
measurements were done on a Quantum Design (QD) 
SQUID MPMS in magnetic fields up to 5 T and using the 
vibrating sample magnetometer (VSM) setup of Quan¬ 
tum Design PPMS up to 14 T in the temperature range 
of 1.6-300 K. Heat capacity data were acquired with the 
QD PPMS in fields up to 14 T. 

Electronic structure calculations within density func¬ 
tional theory (DFT) were performed with the full- 
potential local-orbital code FPL09. 07-41"^° on the 100 K 
crystal structure in combination with the local density 
approximation (LDA),^^ generalized gradient approxi¬ 
mation (GGA )42 DFT+U method.^3.44 ^ 4^4x4 

fc-mesh was employed for LDA and GGA runs while su¬ 
per cells used for DFT-I-U calculations were computed 
for about 100 fc-points in the symmetry-irreducible part 
of the first Brillouin zone. We also performed auxiliary 
calculations using the Heyd-Scuseria-Ernzerhof (HSE06) 
hybrid DFT-functional^® as implemented in the Vienna 
Ab initio Simulation Package (VASP5.2) code.^® 

The hydrogen positions, which are essential for the cal¬ 
culation of the exchange couplings, have not been de¬ 
termined so far^^ since H is almost invisible in XRD due 
to its very low scattering power. Alternative experimen¬ 
tal techniques such as neutron diffraction require large 
and, preferably, deuterated samples that are not avail¬ 
able in nature. Therefore, we determined the positions 
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FIG. 1. (Color online) The left panel shows the crystal structure of langite, Cu 4 ( 0 H) 6 S 04 • 2 H 2 O. CUO 4 plaquettes are shown 
in orange, H are shown as grey spheres and SO 4 groups are shown in yellow. In the right panel, a single crystallographic layer 
is displayed where different colors for the four different Cn-positions are used. The picture in the center shows small light-blue 
langite crystals from the Podlipa and Reinera mine, Lubietova, Slovakia. 


of hydrogen by numerical optimization of the atomic pa¬ 
rameters with respect to a minimization of the total en¬ 
ergy. These calculations were performed within GGA 
and have proved to be highly efficient and sufficiently 
accurate for cuprates in recent studies. 

The exchange coupling constants J^- were calculated 
within DFT following two different strategies. One strat¬ 
egy involves the analysis of the half-filled LDA bands at 
the Fermi level allowing for the determination of leading 
exchange pathways by an evaluation of the electron hop¬ 
ping integrals Uj. The Uj are computed as off-diagonal 
Hamiltonian matrix elements of Cu-centered Wannier 
functions (WFs) constructed for the half-filled bands. 
The spurious metallic state produced for magnetic in¬ 
sulators within LDA can be remedied by inserting LDA- 
based tij into an effective Hubbard model with the ef¬ 
fective onsite Coulomb repulsion UeS, where in cuprates 
typically UeS — 4.5 In the limit of strong cor¬ 

relations, tij <C Heff, which is perfectly fulfilled in langite 
(see Table H), antiferromagnetic (AFM) contributions to 
the total exchange constants Jij can be estimated in sec¬ 
ond order as Jfj^^ = it’^j/UeS- A more detailed descrip¬ 
tion of the procedure can be found, e.g., in Refs. 28 and 
49. 

Alternatively, strong electron correlations are added 
on top of LDA by the LSDA-|-f7 method in a mean-field 
way and are thus included in the self-consistent proce¬ 
dure. This allows for calculating total exchange constants 
which contain also the ferromag¬ 
netic (FM) contributions. A fully localized limit (FLL) 
approximation was used for correcting the double count¬ 
ing. The on-site Coulomb repulsion and Hand’s exchange 
were set to Ud = 8.5± 1 eV and Jd = 1 eV, respectively, a 
choice which has been successfully used for several other 
cuprates.The total exchange coupling constants 


Jij of the spin Hamiltonian 

H = JijSi ■ Sj ( 1 ) 

are calculated as differences between total energies of var¬ 
ious collinear (broken-symmetry) spin states. 

III. CRYSTAL STRUCTURE 

Using the experimental crystal structure of langite re¬ 
ported in Ref. 47, we first routinely performed a DFT- 
optimization of the atomic parameters of all atoms in 
the unit cell with the lattice parameters being fixed 
to their experimental values. Deviations up to 0.3 A 
between the experimental and optimized Cu-0 bond 
lengths prompted us to reinvestigate the crystal struc¬ 
ture of langite with single-crystal XRD. We also per¬ 
formed low-temperature XRD measurements in order to 
probe possible temperature-induced structural changes 
that may be relevant to understanding the magnetism. 

Table I compiles the results of the structural study at 
100 K, which served as input for all DFT calculations. 
Additional crystallographic data collected at 140, 220, 
250, 280 K are provided in the Supplementary Material.®^ 
In the temperature range between 100-280 K the unit cell 
volume increases by about 0.9% with increasing T. The 
largest change in the lattice parameters was observed not 
along a perpendicular to the structural layers (Fig. 1), 
as one might intuitively expect, but along the c direc¬ 
tion. With increasing T, the c parameter increases by 
about 0.35% arising reflecting the flattening of the lay¬ 
ers. Changes along the a and b axes are similar, about 
0.26% each. The monoclinic angle remains almost con¬ 
stant for the investigated temperature range. 

In the presently available structural data,^' hydrogen 
positions have been determined on a semiempirical level. 
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only. One of the hydrogen has been placed on the sul¬ 
phate group which is quite unexpected. In the related 
Cu-sulphate brochantite, Cu 4 ( 0 H)eS 04 , H-atoms have 
been reliably located by neutron diffraction on a deuter- 
ated sample, and no hydrogen was found at the SO 4 
groups but at the Cu-0 layers.^® More doubts on the reli¬ 
ability of the tentative H-positions of langite as provided 
in Ref. 47 arise from the geometry of the water molecules. 
While one of them shows bond lengths close to a free wa¬ 
ter molecule, the other one is strongly distorted with 0-H 
distances of 0.919 A and 1.032 A, respectively, and a H~ 
0-H angle of only 88.54. These issues already call for a 
reinvestigation of the hydrogen positions in langite. Be¬ 
sides gaining new structural information, accurate atomic 
H-positions are also essential for the computation of ex¬ 
change coupling constants which are very sensitive to O- 
H distances and the position of H with respect to the 
Cu 04 plaquette planes (Fig. 1).^^ New atomic hydrogen 
positions are given in Table I, which were obtained by 
GGA-optimization (see Sec. II) using various tentative 
positions as starting values to test the stability of the re¬ 
sults. When only hydrogen atoms were allowed to relax, 
the forces on the oxygen atoms of water molecules (OWl 
and OW2) turned out to be quite large, while one of the 
hydrogen atoms moved towards the SO 4 group. Though, 
such a situation cannot fully be excluded and may arise 
for a certain temperature regime due to the spatial prox¬ 
imity of layers and S04 groups it appears unlikely as 
explained before. In a further step, the positions of all H 
atoms together with those of OWl and OW2 have been 
relaxed. This way, we could stabilize the anticipated lan¬ 
gite structure by 2.6eV/unit cell, while the HSO 4 config¬ 
uration became energetically highly unfavorable. A full 
relaxation of all atomic positions further confirms the sta¬ 
bility of the Cu 4 ( 0 H) 6 S 04 • 2 H 2 O structure and shows 
no signatures of the HSO 4 groups. For the two differ¬ 
ent water molecules, 0-H bond lengths between 0.985- 
0.997A and H-O-H angles of 103.8°and 109.2°, respec¬ 
tively, have been obtained, i.e. there are no asymmetrical 
distortions as proposed in the previous structural work.^^ 
A plot showing the hydrogen bonds and the bonding be¬ 
tween the SO 4 groups, water molecules and Cu-0 lay¬ 
ers can be found in the supplementary material.®^ LDA 
band structures and density of states around the Fermi 
level computed for the different crystal structures, i.e. 
from Ref. 47 and our data collected at 100 K, are pro¬ 
vided in the supplementary material®^. LDA-calculations 
on the crystal structures with optimized hydrogen po¬ 
sitions and optimized positions of oxygen in the water 
molecules (OW) are shown as well. Band shifts between 
50-100 meV and considerable changes in the band dis¬ 
persion are observed particularly between -0.6 to -0.1 eV. 
Since the LDA-bands around the Fermi level crucially 
determine the exchange interactions, these data demon¬ 
strate how crucial hydrogen positions and accurate crys¬ 
tal structures are for computing a microscopic magnetic 
model. 

Table I summarizes atomic positions in langite, includ¬ 


ing the OWl and OW2 positions determined both exper¬ 
imentally and by the GGA-optimization. The difference 
between the experimental and computational positions of 
water molecules may reflect temperature-induced struc¬ 
tural changes, because DFT yields the crystal structure 
at zero temperature, whereas experimental structure de¬ 
termination has been performed down to 100 K, only. 
However, we did not observe any sharp structural phase 
transitions below 100 K in thermodynamic properties re¬ 
ported in Sec. IV. It is also possible that the discrepancy 
between the experimental and computational positions of 
water molecules is intrinsic and related to marginal dis¬ 
order, which is a plausible explanation, given the weak 
(hydrogen) bonding between the water molecules and 
the rest of the crystal structure. Vibration spectroscopy 
could provide further insight into the nature of hydrogen 
bonding and positions of water molecules in langite, but 
it lies beyond the scope of our study, which is focused on 
the magnetism of langite. Relevant magnetic interactions 
run within the Cu-0 layers and should not depend on the 
exact positions of the out-of-plane water molecules. For 
the sake of consistency and given the fact that magnetic 
ordering in langite occurs well below 100 K, we used the 
relaxed positions of OWl and OW2 in the further micro¬ 
scopic analysis (Sec. V). 


The crystal structure of langite features four different 
Cu-positions. The basic building unit are layers formed 
by planar chains of edge-sharing CUO 4 plaquettes (type 
A chain) as well as buckled chains of corner-sharing CUO 4 
plaquettes (type B chains) (Fig. 1), where the chains are 
directly linked to each other. Sulphate groups and water 
molecules are located between the layers. The Gu-0- 
Cu bridging angles, which are of crucial importance for 
the exchange couplings between the Cu-sites, amount to 
99.49°/99.11°, 98.64°/97.81° in the type-A chains. Be¬ 
tween two edgesharing plaquettes in the type-A chains, 
the two bridging angles are different, i.e. the bridge is 
not symmetrical. Both angles are given separated by /•.• 
In the cornersharing type-B chains two different bridging 
angles occur of 101.05° and 104.71° , respectively (see 
Table II) . A figure showing Gu-0 bonding distances and 
bridging angles of the two different chain types is pro¬ 
vided in the supplementary material.®^ The bridging an¬ 
gles between the two chain types, A and B, are between 
105-109°, i.e. the layers are strongly buckled (Fig. 1). 
According to Goodenough-Kanamori rules, one expects 
ferromagnetic (FM) exchange for bridging angles close 
to 90° and antiferromagnetic (AFM) exchange for larger 
bridging angles. The crossover is at about 95-100°,^®’^^ 
and the exchange couplings in the edge-sharing chains 
of langite are difficult to guess in this transition region, 
even qualitatively (see Sec. V), while all other couplings 
would naively be assumed AFM. However, such simple 
considerations are bound to fail for langite as will be 
demonstrated in Secs. V and VI below. 
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TABLE I. Refined atomic positions (in fractions of lattice pa¬ 
rameters) and isotropic atomic displacement parameters Uiso 
(in xlO“^ A^) of langite, Cu 4 ( 0 H) 6 S 04 • 2 H 2 O, collected at 
100 K. Refinement residuals are Ri = 4.64%, WR 2 = 7.6%. 
All atoms are in the general position 2a of the space group 
Plcl. The lattice parameters are as follows: a = 7.1231(8) A, 
b = 6.0305(7) A, c = 11.1935(12) A and /3 = 90.1479(14)°. 
OW and HW denote the O and H atoms of the H 2 O molecules. 
The OS atoms belong to the SO 4 tetrahedra. The H-positions 
have been obtained by numerical optimization within GGA. 
For OW, experimental and optimized positions are provided. 


Atom 

xja 

Vib 

z/c 

U iso 

Gul 

0.99960(13) 

0.99762(18) 

0.49990(9) 

0.46(2) 

Gu2 

0.99260(13) 

0.49213(19) 

0.50188(9) 

0.43(3) 

Gu3 

0.00399(14) 

0.75566(18) 

0.75295(11) 

0.40(2) 

Gu4 

0.00878(15) 

0.25490(17) 

0.75183(11) 

0.38(2) 

S 

0.5778(4) 

0.1854(3) 

0.4201(3) 

0.46(4) 

01 

0.8860(9) 

0.0001(10) 

0.6634(6) 

0.41(12) 

02 

0.8882(9) 

0.5044(10) 

0.6652(6) 

0.54(13) 

03 

0.1156(9) 

0.5098(10) 

0.8432(7) 

0.43(13) 

04 

0.1412(8) 

0.2452(11) 

0.5604(6) 

0.39(13) 

05 

0.8600(8) 

0.7441(9) 

0.4412(6) 

0.45(13) 

06 

0.1256(9) 

0.0054(10) 

0.8398(7) 

0.56(14) 

OSl 

0.7845(8) 

0.2311(9) 

0.4116(5) 

0.42(12) 

OS2 

0.5388(8) 

0.0666(9) 

0.5344(6) 

1.08(13) 

OS3 

0.4770(8) 

0.4015(10) 

0.4219(6) 

0.94(12) 

OS4 

0.5165(8) 

0.9519(10) 

0.8163(6) 

0.91(13) 

OWl 

0.2626(9) 

0.7398(10) 

0.6008(6) 

0.73(13) 

OW2 

0.5178(9) 

0.4283(11) 

0.6955(7) 

1.16(15) 


GGA-optimization 



OWl 

0.27422 

0.73641 

0.59891 


OW2 

0.52514 

0.43403 

0.69849 


HI 

0.27170 

0.74411 

0.03266 


H2 

0.35870 

0.38187 

0.12901 


H3 

0.72674 

0.74906 

0.46285 


H4 

0.26599 

0.00612 

0.33549 


H6 

0.74761 

0.99558 

0.65733 


H7 

0.25162 

0.46853 

0.36133 


HlWl 

0.51463 

0.54313 

0.28622 


H2W1 

0.36074 

0.85755 

0.57697 


H1W2 

0.51392 

0.27508 

0.68068 


H2W2 

0.74644 

0.49444 

0.67040 



IV. THERMODYNAMIC PROPERTIES 

All measurements presented in this section were per¬ 
formed on powder from the same specimen as the one 
used for the single-crystal XRD. The powder quality has 
been diligently checked by powder XRD, revealing almost 
pure langite.®^ 

The temperature-dependent magnetic susceptibility 
x(T) measured in magnetic fields of 0.5 and 2T is shown 
in Fig. 2, where the two curves are almost identical. A 


Curie-Weiss fit x(r) = C/{T+9) of the high-temperature 
regime (220-290 K) of the 2 T data yields 9 = 18.2 K and 
C = 0.481 emu K (mol Cu)“^. From the constant C, 
we obtain an effective magnetic moment of 1.96 ^b/Cu 
which is larger than the spin only value of 1.73 and im¬ 
plies the g-factor of 2.26 lying still in the expected range 
for Cu^+ compounds.The positive 9 indicates pre¬ 
dominant antiferromagnetic couplings, which are, how¬ 
ever, quite weak. In the low-temperature regime, x(’^) 
features a rather sharp peak at 7.5 K. This peak is some¬ 
what asymmetric and thus different from the suscepti¬ 
bility maxima in conventional low-dimensional antifer- 
romagnets, where short-range magnetic order is formed 
well above the Neel temperature While no in¬ 

dications of a magnetic transition are seen in the raw 
susceptibility data, Fisher’s heat capacity d{xT)/dT re¬ 
veals a kink around 5.5 K that can be paralleled to the 
anomaly in the specific heat and ascribed to the mag¬ 
netic ordering transition. The absence of a Curie tail at 
the lowest T, typically arising from paramagnetic spin- 
1/2 impurities (see e.g. Ref. 28), demonstrates the high 
quality of our natural sample. 

The 5.7 K anomaly in the specific heat generally re¬ 
sembles a A-type anomaly, which is expected at a second- 
order phase transition. The broadening of this anomaly 
may be driven by effects of magnetic anisotropy. Mag¬ 
netic nature of the 5.7 K transition is corroborated by 
its field dependence. Despite the relatively low value of 
Tat, the transition is well visible up to at least 14T, and 
the transition temperature changes only slightly in the 
applied magnetic field. At higher temperatures, the lat¬ 
tice contribution to the specific heat dominates. Below 
the transition temperature Tjy, Cmag{T) decreases, but 
it does not follow the simple behavior of a 2D an- 
tiferromagnet which might be anticipated for a layered 
system like langite. However, it is also not following a 
dependence expected for a conventional 3D antiferro- 
magnet (see Fig. 3). 

The magnetic contribution Cmag was obtained by sub¬ 
tracting the lattice contribution Ciat from the measured 
Cp data, where C'iat(T) was approximated by fitting a 

n—7 

polynomial®^ C'iat(T') = X) CnT'^, proposed by Johnston 

n—3 

et to the Cp{T) data in the temperature range of 
20-39 K. The same polynomial was used to extract the 
magnetic contributions from specific heat data measured 
in the various magnetic fields. In zero magnetic field, 
the magnetic entropy S'mag, released within the magnetic 
transition, was estimated to about 6.8 J/(mol K) by inte¬ 
grating Cinag/T. Thus, only about 30% of the expected 
'S'mag = 4i?ln2 for a spin-1/2 system is released within 
the transition anomaly and right above Tjv, while the 
rest is spread towards higher temperatures, which is typ¬ 
ical for low-dimensional antiferromagnets®® and corrob¬ 
orates that Tat is somewhat lower than the energy scale 
of the exchange couplings given by, e.g., 0 ~ 18K. A 
similar value for S'mag has been reported for the related 
mineral brochantite (see also Sec. VI) releasing 7.9 J/(mol 
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FIG. 2. (Color online) The experimental susceptibility data 
x{T) of langite, Cu 4 (OH) 6 S 04 • 2 H 2 O, collected at magnetic 
fields of 0.5 and 2.0 T on a powder sample. d{xT)/dT is shown 
in the inset as blue line. 

K), which is about 34% of the total magnetic entropy, in 
the vicinity of the magnetic transition.^® 

Field-dependent magnetization M(H) (Fig. 4) mea¬ 
sured in fields up to 14 T features a kink around 4T, 
reaches half-saturation around 12 T, and keeps increas¬ 
ing up to at least 14 T. The kink at 4T is reminiscent of 
a spin-flop transition that, however, happens at a much 
higher field than in other Cu^'*' oxides.®^ Above 2 K, the 
features of the magnetization curve are smeared out, so 
we were not able to map them as a function of tempera¬ 
ture and construct a comprehensive T—H phase diagram. 

V. ELECTRONIC STRUCTURE AND 

MAGNETIC EXCHANGE COUPLINGS 

In this section, we derive a microscopic magnetic model 
that could be used to understand the complex behavior 
of langite. Microscopic models based on empirical con¬ 
siderations are prone to error because superexchange in 
Cu^'*' compounds depends on tiny structural details and 
cannot be fully captured by empirical rules. Moreover, 
the presence of four distinct Cu sites in the crystal struc¬ 
ture implies that interactions with similar Cu-Cu dis¬ 
tances and superexchange pathways are not related by 
symmetry and may be unequal. Therefore, an empirical 
approach for deriving a microscopic magnetic model is 
bound to fail for langite. Accordingly, we employ numer¬ 
ical electronic structure DFT calculations allowing for 
a direct computation of individual exchange couplings 
Jij. In combination with numerical simulations of the 
thermodynamical properties, such calculations often pro¬ 
vided consistent description of the macroscopic magnetic 
behavior based on microscopic considerations.^®’^®’®®’®® 

With suitably chosen correlation parameters, such as 
the Coulomb repulsion Ud in LSDA-bU, one expects that 



FIG. 3. (Golor online) The magnetic contribution to the spe¬ 
cific heat, Cmag(T), of langite, Gu 4 ( 0 H) 6 S 04 • 2 H 2 O, mea¬ 
sured in magnetic helds between 0-14 T on a powder sample. 
The inset shows a double logarithmic plot of the magnetic 
specific heat Cmag{T) of langite collected in zero magnetic 
field. The green and blue lines show T® and T® fits, respec¬ 
tively. T® and T® would be the anticipated behaviors of a 
two-dimensional and three-dimensional antiferromagnet, re¬ 
spectively. The grey dotted curve shows the fitted lattice 
background Ciat(T). 



FIG. 4. (Golor online) Field-dependent magnetization data 
M(H) of langite, Cu 4 ( 0 H) 6 S 04 • 2 H 2 O, up to 14T collected 
at a temperature of 2.OK on a powder sample. The arrows 
indicate the kink at about 4T and half saturation which is 
reached at about 12 T assuming a g-factor of 2.0. 


DFT results are accurate within 10% for most insulating 
spin-1/2 materials and the respective interaction path¬ 
ways. However, the error bars increase for very weak 
couplings and for those couplings, where special nature 
of the superexchange pathway renders ferro- and anti¬ 
ferromagnetic contributions comparable in size. Further 
information on the computational procedure and the ac- 
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curacy of computed exchange couplings can be found in 
Refs. 13, 29, 59-62. 

As a first step, LDA calculations were performed, 
yielding a broad valence band complex of about 10 eV 
(Fig. 5), which is typical for cuprates.Low-energy 
magnetic excitations should be largely determined by the 
band complex of eight half-filled bands around the Fermi 
level, between —0.5 and 0.45 eV. The eight bands arise 
from the eight Cu^“'"-ions per unit cell and their corre¬ 
sponding eight half-hlled Sd-orbitals. Local coordinate 
systems on the eight Cu-sites (with the local z-axis cho¬ 
sen perpendicular to the Cu 04 -planes and the local x- 
axis oriented parallel to a Cu-0 bond) allow analyzing 
the orbital character of the half-filled bands. They are es¬ 
sentially of Cu( 3 (ia; 2 _j, 2 ) and 0{2px,2py) character while 
admixtures from H 2 O and particularly SO 4 groups are 
small. Accordingly, the latter molecules do not play a di¬ 
rect role for the exchange couplings in langite. This set of 
eight bands is now projected onto Cu-centered Wannier 
functions (WFs) to evaluate the hopping parameters tij. 
Owing to the four different Cu-sites, many different ex¬ 
change pathways are effective in langite. Table II lists all 
\tij\ > 20meV as well as the corresponding The 

largest interlayer hopping f is only about — 5meV ren¬ 
dering the spin lattice of langite nearly two-dimensional. 
The positions of the respective exchange pathways in the 
crystal structure are shown in Fig. 6 . 

Full exchange constants Jij = J™ -I- computed 

with the LSDA-I-C/ method are provided in the last col¬ 
umn of Table II, where error bars show the effect of 
changing the Coulomb repulsion parameter by ±1 eV. 
This parameter affects absolute values of exchange cou¬ 
plings, while their ratios typically change by few percent 
only. However, for weak couplings error bars can exceed 
absolute values, and thus the ratios are strongly affected 
as well. 

The two nearest-neighbor (NN) couplings of the Cu- 
spins on the type-B chains, Js and J(, are both FM, 
with being much smaller than Jg. This results in fer¬ 
romagnetic spin chains with alternating exchange cou¬ 
plings. These spin chains interact antiferromagnetically 
via Jjj and represent the first magnetic sublattice (SLl). 
The difference between Jg and J( can be traced back to 
the relevant Cu-O-Cu angles (Table II) . The smaller an¬ 
gle for Jg leads to a stronger FM interaction (see also 
Sec. VI). 

The NN coupling Ji of the Cu-spins on the edge¬ 
sharing type-A chains is AFM and about four times 
stronger than the weak coupling J[. The AFM next- 
nearest-neighbor (NNN) couplings J 2 and J 2 are of the 
same strength as Ji. The exchange interactions of the 
spins on the type-A chains may, thus, be described in 
terms of magnetic two-leg ladders (Fig. 6 ), where J 2 and 
J 2 build the legs, while Ji forms the rungs, and J( is 
a frustrating diagonal interaction. This represents the 
second magnetic sublattice (SL2). The strengths of Ji 
and J[ can be again traced back to the Cu-O-Cu angles. 
The smaller bridging angles render J( weaker than Ji. 



FIG. 5. (Color online) The top panel shows the total and 
partial density of states (DOS) from LDA calculations. In 
the lower panel the eight half-filled LDA bands around the 
Fermi level are shown. ’’Wannier” denotes bands calcnlated 
with Cu-centered Wannier functions. 

However, these antiferromagnetic couplings are observed 
for the bridging angles below 100 °, while the ferromag¬ 
netic couplings Js and J' occur for the bridging angles 
above 100°. This instructive situation highlights limita¬ 
tions of the Goodenough-Kanamori-Anderson rules and 
the importance of the mutual arrangement of the CUO 4 
plaquettes, which share edges (Ji and J[) or corners {Js 
and J(), respectively. 

Multiple couplings between the two sublattices are 
mostly weak. The strongest inter-sublattice interaction 
Jh is comparable in size to Ji, J 2 , and J 2 . However, it 
is less abundant than the intra-sublattice couplings, and 
on average one finds that only half of Jh contributes to 
the effective molecular held on the SL2. Therefore, in 
a hrst approximation one can consider langite as anti- 
ferromagnet built of two sublattices, where the sublat¬ 
tice SL2 is ID, the sublattice SLl is 2D, and the inter¬ 
sublattice couplings are weaker than the leading cou¬ 
plings within each of the sublattices. The interlayer cou¬ 
pling J', which is about 0.2K, may be responsible for 
the long-range magnetic order observed in langite be¬ 
cause three-dimensional order requires the coupling be¬ 
tween the layers (Fig. 3). We refrained from estimating 
J' using LSDA-I-C/, though, because such small couplings 














































FIG. 6 . (Color online) The left panel shows the intralayer exchange pathways. The central panel shows the two structural 
chain types and the respective intrachain exchange pathways. The exchange interactions of Cu-spins on the type-B chains may 
be described in terms of alternating ferromagnetic spin chains which are coupled with each other by J//. This represents the 
first magnetic sublattice (SLl). A two-leg spin ladder drawn from the intrachain exchange couplings of the Cu-spins on the 
type-A chains is shown in the right panel. These antiferromagnetic spin ladders form the second magnetic sublattice (SL2). A 
list of all exchange couplings, their strengths and structural characteristics are provided in Table II. 


are hard to calculate reliably.'*’® 

The exchange couplings given in Table II allow estimat¬ 
ing a ” local Curie-Weiss temperature” 9k for each Cu-site 
i according to 9k = S{S -|- l)/3 • ' Ji) where S is the 

i 

electron spin and Zi shows how often a certain coupling Ji 
occurs on a given site. 9k is thereby a measure for the lo¬ 
cal coupling strengths on the Cu-site k. Accordingly, we 
get for the four Cu-sites (Fig. 6 ): 9\ = — 15K, 02 = —3K, 
03 = 31K, 04 = 27 K. The overall Curie-Weiss temper¬ 
ature may be approximated by averaging over all sites 
yielding 0 = 10 K which is in reasonable agreement with 
the experimental value of 18.2 K (see Sec. IV) regarding 
the intricate microscopic magnetic model and the large 
number of exchange couplings. The striking difference 
between 0 i and 02 arises from the nearest-neighbor in¬ 
terchain couplings, where the ferromagnetic Ja operates 
on Cu-site 1 and the antiferromagnetic Jh is effective on 
Cu-site 2. It is worth noting that the single experimental 
parameter 0 is usually sufficient for verifying the choice of 
Ud in LSDA-bC/, because the uncertainty in absolute val¬ 
ues of the computed exchange couplings is much higher 
than in their ratios, and thus only the absolute scale of 
computed exchange couplings should be cross-checked ex¬ 
perimentally. On the other hand, experimental evalua¬ 
tion of individual exchange couplings in langite may be 
an arduous task, given the overall complexity of the spin 
lattice. In Sec. VI, we further discuss ramifications of our 
microscopic magnetic model and its relevant macroscopic 
features that can be tracked experimentally. 

We also performed calculations using a HSE06 hybrid 
DFT-functional,^® as implemented in VASP5. 2.^® These 


calculations were feasible for short-range couplings only. 
The long-range couplings would require big supercells 
that cannot be treated with the computationally expen¬ 
sive HSE06 method with the required accuracy. 

In contrast to LSDA-|-[/, the hybrid-functional ap¬ 
proach does not include the effect of local Coulomb 
repulsion explicitly, thus leading to less accurate esti¬ 
mates of individual exchange couplings and the overesti¬ 
mate of ferromagnetic terms.'*® On the other hand, this 
method is free from adjustable parameters and does not 
involve the ambiguous choice of the Coulomb repulsion 
Ud- The HSE06 results can be found in the supplemen¬ 
tary material.®* They are generally similar to those from 
LSDA-|-f7 (Table II) and confirm main features of the 
langite spin lattice: i) the FM nature of Js and ii) 
the AFM nature of Ji and J[ , and iii) the | Jg | > | J' | and 
J[ < Ji trends discussed above. The inter-sublattice cou¬ 
plings are mostly FM in HSE06 because of the general 
tendency of hybrid functionals to overestimate ferromag¬ 
netic contributions to the exchange. 

Disregarding the weak coupling J[, we find that indi¬ 
vidual magnetic sublattices of langite are non-frustrated 
because none of the leading couplings Jg, J', and J// 
for the type-A chains and Ji, J 2 , and J 2 for the type-B 
chains compete to each other. The two leading couplings 
between the sublattices, FM Ja and AFM Jh, are not 
frustrated either, because they are compatible with the 
AFM order between the FM type-B chains, as imposed by 
J 77 . However, other inter-sublattice couplings frustrate 
the spin lattice and render it very complex. Given the 
large number of non-equivalent exchange couplings and 
their frustrated nature, we restrict ourselves to a qualita- 
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TABLE II. The transfer integrals tij (meV) and the AFM 
contributions to the exchange constants = AtljlUeS 

(K) where ?7efi = 4.5eV. d(Cu-Cu) and Cu-O-Cu denote Cu- 
Cu distances (A) and Cu-O-Cu angles (deg), respectively. 
The Jij (K) given in the last column are calculated with the 
LSDA+17 method and Ud = 8.5±1 eV and Jd = 1 eV. The dif¬ 
ferent groups of exchange couplings are sorted with respect to 
Cu-Cu distances. For a detailed explanation of the bridging 
angles see Sec. Ill as well as the supplementary material.®^ 



d(Cu-Cu) Cu-O-Cu 

tij 

tAFM 

•Jij 

Jij 



type A chains (edge-sharing) 


J'l 

3.011 

97.81/98.64 

-144 

213 

9± 15 

Jl 

3.020 

99.10/99.49 

-155 

247 

38 ±20 

J2 

6.030 


61 

39 

36 ±6 

J2 

6.030 


59 

36 

35 ±6 



type B chains (corner-sharing) 


Js 

2.983 

101.05 

40 

16 

-74 ± 10 

J's 

3.049 

104.706 

62 

40 

-23 ±3 

Jb 2 

6.030 


31 

10 

< 1 

J's 2 

6.030 


30 

9 

< 1 


nearest-neighbor interchain couplings 


Ja 

3.139 

103.55 

85 

74 

-19 ± 1 

■h 

3.144 

105.08 

-109 

123 

5±5 

Jc 

3.163 

106.64 

91 

85 

-6 ±2 

Jd 

3.168 

105.23 

109 

121 

8±5 

Je 

3.186 

106.65 

-91 

85 

0±4 

Ji 

3.190 

106.61 

94 

91 

-12 ±2 

J9 

3.219 

106.10 

-92 

87 

-13 ± 1 

Jh 

3.229 

109.23 

-107 

119 

31 ±7 


next-nearest-neighbor interchain 

couplings 


Jl 

6.309 


-51 

27 

10 ± 10 

Jll 

6.407 


-61 

38 

51 ±20 



interlayer coupling 



J' 

7.874 


-5 

0.2 



five discussion of the magnetic behavior in Sec. VI below. 


VI. DISCUSSION 

Magnetic properties of the natural Cu^“''-mineral lan- 
gite are peculiar, yet complicated. Its crystal structure 
consists of layers formed by directly connected and alter¬ 
nating ordered edge- and corner-sharing chains of CUO 4 
plaquettes. Such structural motives can be found in sev¬ 
eral compounds (see Sec. Ill) that have same topology 
of the magnetic layer, but slightly different Cu-O-Cu 
angles and, thus, potentially different exchange scenar¬ 
ios. Magnetic properties of only one of these materials, 
brochantite, have been reported. Therefore, it is still 
an open question which magnetic properties arise when 
edge- and corner-sharing Cu^“'" chains are joined into lay¬ 
ers, and how these properties are affected by structural 


details. 

Our study shows that such layers cannot be viewed 
as a simple combination of weakly coupled magnetic 
chains. Unanticipated interchain couplings, such as J// 
that features an unusually long superexchange pathway, 
render the spin lattice much more complex. Remarkably, 
though, we can still split this lattice into two sublat¬ 
tices composed of type-A and type-B chains, respectively. 
The difference between these sublattices and individual 
chains pertains to the fact that all type-B chains form 
a single 2D sublattice SLI, whereas sublattice SL2 com¬ 
prises weakly coupled type-A chains and thus remains 
effectively ID. We should also emphasize that the inter¬ 
sublattice couplings are clearly non-negligible. A quan¬ 
titative description of langite will, therefore, require the 
consideration of the full spin lattice that is partially frus¬ 
trated. This problem must be tackled with advanced 
simulation techniques and lies beyond the scope of our 
present study, where we restrict ourselves to the quali¬ 
tative analysis and demonstrate that the model of two 
different sublattices can be used to rationalize main fea¬ 
tures of the experimental data. 

First, the abrupt increase of the magnetization in low 
fields and the fact that half-saturation is reached already 
at 10 — 12 T is consistent with the presence of sublat¬ 
tice SLI, which is largely ferromagnetic. This sublattice 
should saturate as soon as magnetic field overcomes the 
effect of the AFM coupling Jjj. There is only one J// 
coupling per Cu site, so the half-saturation should be 
reached at Hgi = kBJii/— 38T, which is much 
higher than 12 T observed experimentally. The origin 
of this discrepancy is not entirely clear. The presence 
of the second magnetic sublattice SL2 (comprising the 
spins on the type-A chains) and the frustrated interac¬ 
tions between the sublattices may overcome the effect of 
J 77 and facilitate the saturation of the SLI (consisting 
of spins on the type-B chain) already in low field, al¬ 
though a detailed investigation of this behavior requires 
numerical simulations for the full spin lattice of langite, 
which are not feasible, as we explained above. The SL2 
is antiferromagnetic and its saturation is expected at 
Hs 2 = kB/{ggB){Ji + J 2 + J 2 ) — 81T which would 
be interesting to probe experimentally. We thus expect 
that above 14 T the magnetization of langite increases 
much slower than in low fields, and that full saturation 
is reached around 80 T reflecting the presence of sizable 
AFM couplings in this system. 

Magnetic susceptibility of langite lacks a broad maxi¬ 
mum that would be expected in a quasi-2D antiferromag- 
net. This observation is also consistent with the presence 
of the mostly FM sublattice SLI (Fig. 6 ) which lacks any 
susceptibility maximum down to and, thus, masks the 
susceptibility maximum related to the AFM sublattice 
SL2. A similar behavior has been observed in CuP 206 ,^® 
where none of the sublattices is ferromagnetic, but very 
weak couplings in one of the sublattices render half of 
the spins paramagnetic down to low temperatures, and 
no susceptibility maximum is observed down to Tjv. It is 
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worth noting that the asymmetric maximum in the sus¬ 
ceptibility of langite around 8 K cannot be taken as a typ¬ 
ical signature of short-range order in a quasi-2D system, 
because this maximum is observed at temperatures well 
below the Curie-Weiss temperature 9 ~ 18 K, whereas 
in a 2D system, e.g., in a square-lattice antiferromagnet, 
T'nax ~ 0 is expected. 

The temperature of the antiferromagnetic ordering in 
langite, T^v — 5.7K, is quite low compared to leading ex¬ 
change couplings I Js|, Ji, Jii, J 2 , and J 2 that are at least 
35 — 40 K each. We tentatively find J < 0.2, which 
is very low for a quasi-2D antiferromagnet,®^ although 
a correct definition of an effective intralayer coupling J 
may be difficult in this case given the very complex na¬ 
ture of the spin lattice. 

A frustration ratio of O/Tj^ ~ 3 is less impressive, but 
one has to acknowledge that the macroscopic d is a sum 
of FM and AFM couplings (Table II) and thus underes¬ 
timates the overall energy scale of exchange couplings in 
langite. 

The magnetic ground state of langite may be peculiar. 
In sublattice SLl, one expects FM order along b and 
AFM order along a and c, arising from the interchain 
interaction J 77 and the weak interlayer coupling J' . The 
sublattice SL2 is a two-leg spin ladder and, when taken on 
its own, features a spin-singlet ground state without long- 
range magnetic order. Although interchain couplings and 
the couplings to SLl will trigger the formation of ordered 
moments even in SL2, these moments are expected to be 
much smaller than in SLl. This difference in the ordered 
magnetic moments is one of the fingerprints of the two- 
sublattice model and can be probed experimentally by 
nuclear magnetic resonance or neutron scattering. 

The drastic difference between the ordered magnetic 
moments on different Cu-sites has been previously seen 
in other Cu^'*' minerals. In antlerite Cu 3 ( 0 H) 4 S 04 ,®^’®^ 
two side chains of the B-type encompass the central chain 
of the A-type that together form a ribbon, which is some¬ 
times considered as a three-leg spin ladder. Neutron 
scattering revealed ordered magnetic moment of 0.88 
on the terminal (type-B) chains and zero magnetic mo¬ 
ment on the central (type-A) chains.®^ A similar type 
of magnetic order is expected in langite, where spins 
in the type-B chains will form long-range magnetic or¬ 
der, whereas spins in the type-A chains should develop a 
gapped ground state with zero ordered moment, as typ¬ 
ical for two-leg spin ladders. This unusual, partial mag¬ 
netic order may be reflected in magnetic excitations and 
macroscopic properties such as specific heat below T^r. 
Indeed, the specific heat of langite clearly deviates from 
the standard or behaviors and remains an inter¬ 
esting problem for future investigation. 

Another Cu^+ mineral, brochantite Cu 4 ( 0 H) 6 S 04 , is 
remarkably different from both langite and antlerite. 
From the chemistry perspective, it is a dehydrated ver¬ 
sion of langite featuring same type of magnetic layers. 
However, details of their geometry are somewhat differ¬ 
ent because water molecules are missing, and the sep¬ 


aration between the layers is about twice shorter than 
in langite. Neutron diffraction reports very small mag¬ 
netic moments within the corner-sharing type-B chains 
( 0 . 22 /is) and much larger ordered moments within the 
type-A chains ( 0 . 74 / 7 ^).®® This is very different from 
the ground state of antlerite (and, presumably, of lan¬ 
gite) and may indicate a different exchange topology. In¬ 
deed, the Curie-Weiss temperature of brochantite {0 ~ 
90 K)®® is much higher than 18 K and 4K in langite and 
antlerite,®® respectively. Moreover, brochantite features 
a broad susceptibility maximum around 60 K, far above 
Tjv, while neither langite nor antlerite show such broad 
maxima. These features suggest that magnetic interac¬ 
tions in brochantite are predominantly AFM, whereas 
langite and antlerite reveal a subtle interplay of FM and 
AFM exchange couplings. Further microscopic insight 
into these differences is clearly needed and requires a sys¬ 
tematic computational study of the aforementioned Cu^'*' 
minerals. 

Naively, the AFM nature of brochantite can be as¬ 
cribed to the larger Cu-O-Cu angles in the range 105.6 — 
122.5°, while the Cu-O-Cu angles in langite are, gener¬ 
ally, smaller (103.5 — 109.2°). However, this simple analy¬ 
sis in the spirit of Goodenough-Kanamori-Anderson rules 
may often be misleading. 

We have mentioned that the values of the Cu-O-Cu 
bridging angles account for \Jg\ > |Js'| and Ji > J(, 
but they do not explain why Jg and J( are FM, while 
Ji and J[ with the smaller bridging angles are AFM. 
Other effects are obviously important in this case. In 
particular, hydrogen atoms bonding to the bridging oxy¬ 
gen have strong influence on the superexchange.Hy¬ 
drogen atoms located out of the CUO 4 planes, favor FM 
exchange. This is definitely relevant for Jg and J( in 
langite where the 0-H bond on the bridging oxygen and 
the Cu 04 plaquettes enclose angles up to 60°. For the 
Cu-mineral clinoclase, Cu 3 (As 04 )( 0 H )3 it was recently 
demonstrated^^ that such a large out-of-plane angle can 
drive the exchange coupling even from a strongly AFM to 
the FM coupling regime. Quite similar results for the in¬ 
trachain physics of the type-B chains have been reported 
in recent studies on antlerite, Cu 3 ( 0 H) 4 S 04 .®®’®^ The 
crystal structure of this compound features triple chains 
consisting of a central type-A chain and type-B chains 
bonded to it on each side. For the Cu-spins on the type- 
B chains an alternating FM coupling has been reported 
from neutron experiments with an antiparallel order be¬ 
tween the chains.®®’®'*. For brochantite, Cu 4 ( 0 H) 6 S 04 , 
featuring structural layers similar to those in langite, neu¬ 
tron data also revealed a FM coupling of the Cu-spins 
within the type-B chains.®® Eventually, in a joint exper¬ 
imental and theoretical study®^ on the rare Cu-mineral 
szenicsite, Cu 3 Mo 04 (OH) 4 , an alternating FM coupling 
on the type-B chains has been reported. This compound 
features triple chains similar to those in antlerite. 
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VII. SUMMARY 

In summary, structural and magnetic properties of the 
Cu^'*' mineral langite have been investigated in a joint 
experimental and theoretical study. Crystal structure of 
langite was refined in the 100 — 280 K temperature range 
using single-crystal XRD, and the H-positions were sub¬ 
sequently determined for the 100 K structure using the¬ 
oretical DFT-approach. The crystal structure consists 
of two types of directly connected Cu-chains, edge- and 
corner-sharing, which form layers separated from each 
other by about 7.5 A. These layers are a common struc¬ 
tural motive in cuprate minerals, but their relevant mag¬ 
netic interactions and resulting magnetic properties have 
been only scarcely investigated. Along with the fact that 
both chain-types taken on their own have revealed fasci¬ 
nating magnetic properties, it intrigued us what kind of 
physics may arise from their combination into layers. 

Our density-functional calculations show that such lay¬ 
ers can not be viewed as a stack of weakly coupled mag¬ 
netic chains. While different chains form different mag¬ 
netic sublattices, interactions between the chains are non- 
negligible, and even the two-sublattice model describes 
the magnetic behavior only qualitatively. It does, how¬ 
ever, capture the crucial feature that the sublattice B 
is predominantly ferromagnetic and prone to the forma¬ 
tion of the long-range order, whereas the sublattice A 
is entirely antiferromagnetic and gapped because of its 
two-leg-ladder geometry. Therefore, we expect a peculiar 
magnetic ground state with drastically different ordered 
moments in the two sublattices. This ground state can be 
paralleled to that of antlerite, where the “idle-spin” be¬ 
havior (no detectable ordered moment) on type-A chains 
has been observed. 

Experimentally, langite undergoes long-range mag¬ 
netic ordering, but at the Neel temperature Tn ~ 5.7 K 
that is well below the Curie-Weiss temperature 0 ~ 18K. 
An effective “frustration ratio” 0/Tn ~ 3 demonstrates 
that the magnetic order in langite is impeded. However, 
the Curie-Weiss temperature is a sum of ferromagnetic 


and antiferromagnetic couplings and thus underestimates 
the energy scale of magnetic exchange. Taking computed 
J’s from Table II, one finds that the Neel temperature of 
langite is remarkably low for a quasi-2D antiferromagnet. 
This reduced value of Tjv is a signature of strong quan¬ 
tum fluctuations that have three concurrent origins: i) 
spin-i nature of Cu^"*" and magnetic low-dimensionality; 
ii) in-plane frustration; iii) proximity of sublattice A to 
the spin-singlet state without long-range magnetic order. 

Altogether, langite is a frustrated quasi-2D antiferro¬ 
magnet that reveals interesting manifestations of quan¬ 
tum magnetism and a peculiar two-sublattice structure 
of the spin lattice. Its ground state is of particular in¬ 
terest for future studies, given the anticipated difference 
between the ordering processes in two magnetic sublat¬ 
tices. Specific heat of langite measured in the ordered 
state does not follow conventional or behavior, 
thus providing first evidence for the unconventional na¬ 
ture of the magnetic ground state and calling for further 
investigation of this interesting material. 


ACKNOWLEDGMENTS 

We are grateful to S.-L. Drechsler and S. Nishimoto 
for valuable discussions. We acknowledge the experimen¬ 
tal support by Yurii Prots, Horst Borrmann (laboratory 
XRD) and Christoph Klausnitzer (low-temperature spe¬ 
cific heat). We would also like to thank Prof. Klaus Thal- 
heim and the Senckenberg Naturhistorische Sammlung 
Dresden for providing the langite sample from Slovakia 
(inventory number 18940 Sy). OJ was partly supported 
by the European Research Council under the European 
Unions Seventh Framework Program FP7/ERC through 
grant agreement n.306447. AT was supported by the Mo- 
bilitas program of the ESF (grant number MTT77) and 
by the Alexander von Humboldt Foundation and Federal 
Ministry for Education and Research, Germany, through 
the Sofja Kovalevskaya Award. 


* st.lebernegg@gmail.com 
^ Helge.Rosner@cpfs.mpg.de 

^ M. Base, I. Terasaki, and K. Uchinokura, Phys. Rev. Lett. 
70 , 3651 (1993). 

^ S. Seki, X. Z. Yu, S. Ishiwata, and Y. Tokura, Science 13 , 
198 (2012). 

® T. Giamarchi, C. Riiegg, and O. Tchernyshyov, Nat. Phys. 
4, 198 (2008). 

J. Schlappa, K. Wohlfeld, K. J. Zhou, M. Mourigal, M. W. 
Haverkort, V. N. Strocov, L. Hozoi, C. Monney, S. Nishi¬ 
moto, S. Singh, A. Revcolevschi, J.-S. Caux, L. Patthey, 
H. M. Ronnow, J. van den Brink, and T. Schmitt, Nature 
485, 82 (2012). 

® H. Rosner, H. Eschrig, R. Hayn, S.-L. Drechsler, and 
J. Malek, Phys. Rev. B 56, 3402 (1997). 


® A. Moller, M. Schmitt, W. Schnelle, T. Forster, and 
H. Rosner, Phys. Rev. B 80, 125106 (2009). 

^ K. Nawa, T. Yajima, Y. Okamoto, and Z. Hiroi, Inorg. 
Chem 54, 5566 (2015). 

® T. Masuda, A. Zheludev, A. Bush, M. Markina, and 

A. Vasiliev, Phys. Rev. Lett. 92, 177201 (2004). 

® L. Capogna, M. Reehuis, A. Maljuk, R. K. Kremer, 

B. Ouladdiaf, M. Jansen, and B. Keimer, Phys. Rev. B 
82, 014407 (2010). 

F. Sapiiia, J. R. guez Carvajal, M. J. Sanchis, R. I. nez, 
A. Beltran, and D. Beltran, Solid State Comm. 74 , 779 
(1990). 

K. Caslin, R. K. Kremer, F. S. Razavi, A. Schulz, 
A. Munoz, F. Pertlik, J. Liu, M.-H. Whangbo, and J. M. 
Law, Phys. Rev. B 89, 014412 (2014). 



12 


M. G. Banks, R. K. Kremer, C. Hoch, A. Simon, B. Oulad- 
diaf, J.-M. Broto, H. Rakoto, C. Lee, and M.-H. Whangbo, 
Phys. Rev. B 80, 024404 (2009). 

C. Lee, J. Liu, M.-H. Whangbo, H.-J. Koo, R. K. Kremer, 
and A. Simon, Phys. Rev. B 86, 060407 (2012). 

B. Willenberg, M. Schapers, K. C. Rule, S. Siillow, 
M. Reehuis, H. Ryll, B. Klemke, K. Kiefer, W. Schotten- 
hamel, B. Buchner, B. Ouladdiaf, M. Uhlarz, R. Beyer, 

J. Wosnitza, and A. U. B. Wolter, Phys. Rev. Lett. 108, 
117202 (2012). 

S. Seki, Y. Yamasaki, M. Soda, M. Matsuura, K. Hirota, 
and Y. Tokura, Phys. Rev. Lett. 100, 127201 (2008). 

Y. Yasui, K. Sato, Y. Kobayashi, and M. Sato, J. Phys. 
Soc. Jpn. 78, 084720 (2009); Y. Kobayashi, K. Sato, Y. Ya¬ 
sui, T. Moyoshi, M. Sato, and K. Kakurai, ibid. 78, 084721 
(2009). 

Y. Yasui, M. Sato, and 1. Terasaki, J. Phys. Soc. Jpn. 80, 
033707 (2011). 

M. Schapers, A. U. B. Wolter, S.-L. Drechsler, S. Nishi- 
moto, K.-H. Muller, M. Abdel-Hafiez, W. Schottenhamel, 
B. Buchner, J. Richter, B. Ouladdiaf, M. Uhlarz, R. Beyer, 
Y. Skourski, J. Wosnitza, K. C. Rule, H. Ryll, B. Klemke, 

K. Kiefer, M. Reehuis, B. Willenberg, and S. Siillow, Phys. 
Rev. B 88, 184410 (2013). 

M. Enderle, B. Fak, H.-J. Mikeska, R. K. Kremer, 
A. Prokofiev, and W. Assmus, Phys. Rev. Lett. 104, 
237207 (2010); S.-L. Drechsler, S. Nishimoto, R. O. 

Kuzian, J. Malek, W. E. A. Lorenz, J. Richter, J. van den 
Brink, M. Schmitt, and H. Rosner, Phys. Rev. Lett. 106, 
219701 (2011); M. Enderle, B. Fak, H.-J. Mikeska, and 
R. Kremer, Phys. Rev. Lett. 106, 219702 (2011). 

M. Mourigal, M. Enderle, B. Fak, R. K. Kremer, J. M. 
Law, A. Schneidewind, A. Hiess, and A. Prokofiev, Phys. 
Rev. Lett. 109, 027203 (2012). 

F. C. Chou, A. Aharony, R. J. Birgeneau, O. Entin- 
Wohlman, M. Greven, A. B. Harris, M. A. Kastner, Y. J. 
Kim, D. S. Kleinberg, Y. S. Lee, and Q. Zhu, Phys. Rev. 
Lett. 78, 535 (1997). 

Y. J. Kim, A. Aharony, R. J. Birgeneau, F. C. Chou, 
O. Entin-Wohlman, R. W. Erwin, M. Greven, A. B. Harris, 
M. A. Kastner, 1. Y. Korenblit, Y. S. Lee, and G. Shirane, 
Phys. Rev. Lett. 83, 852 (1999). 

M. A. Kastner, A. Aharony, R. J. Birgeneau, F. C. Chou, 
O. Entin-Wohlman, M. Greven, A. B. Harris, Y. J. Kim, 
Y. S. Lee, M. E. Parks, and Q. Zhu, Phys. Rev. B 59, 
14702 (1999). 

Y. J. Kim, R. J. Birgeneau, F. C. Chou, M. Greven, M. A. 
Kastner, Y. S. Lee, B. O. Wells, A. Aharony, O. Entin- 
Wohlman, 1. Y. Korenblit, A. B. Harris, R. W. Erwin, 
and G. Shirane, Phys. Rev. B 64, 024435 (2001). 

A. B. Harris, A. Aharony, O. Entin-Wohlman, 1. Y. Ko¬ 
renblit, R. J. Birgeneau, and Y.-J. Kim, Phys. Rev. B 64, 
024436 (2001). 

R. Nath, K. M. Ranjith, J. Sichelschmidt, M. Baenitz, 
Y. Skourski, F. Alet, 1. Rousochatzakis, and A. A. Tsirlin, 
Phys. Rev. B 89, 014407 (2014). 

S. Lebernegg, A. A. Tsirlin, O. Janson, and H. Rosner, 
Phys. Rev. B 87, 235117 (2013). 

S. Lebernegg, A. A. Tsirlin, O. Janson, and H. Rosner, 
Phys. Rev. B 88, 224406 (2013). 

H. Jeschke, 1. Opahle, H. Kandpal, R. Valenti, H. Das, 

T. Saha-Dasgupta, O. Janson, H. Rosner, A. Briihl, 

B. Wolf, M. Lang, J. Richter, S. Hu, X. Wang, R. Pe¬ 
ters, T. Pruschke, and A. Honecker, Phys. Rev. Lett. 106, 


217201 (2011). 

M. Schmitt, O. Janson, M. Schmidt, S. Hoffmann, 

W. Schnelle, S.-L. Drechsler, and H. Rosner, Phys. Rev. 
B 79, 245119 (2009). 

S. Vilminot, M. Richard-Plouet, G. Andre, D. Swierczyn- 
ski, M. Guillot, F. Bouree-Vigneron, and M. Drillon, J. 
Sol. Stat. Chem. 170, 255 (2003). 

S. Lebernegg, A. A. Tsirlin, O. Janson, 1. Rousochatzakis, 
and H. Rosner, unpublished. 

D. Y. Pushcharovsky, R. K. Rastsvetaeva, and H. Sarp, 
J. Alloy. Compd. 239, 23 (1996). 

D. Ginderow and F. Cesbron, Acta Cryst. C 39, 1605 
(1983). 

G. Giester, B. Rieck, and F. Brandstaetter, Miner. Petrol. 
63, 19 (1998). 

®® S. Vilminot, M. Richard-Plouet, G. Andre, D. Swierczyn- 
ski, F. Bouree-Vigneron, and M. Kurmoo, Dalton Trans. 
, 1455 (2006). 

M. R. Bissengaliyeva, N. S. Bekturganov, D. B. Gogol, 
S. T. Taimassova, T. A. Koketai, and M. A. Bespyatov, 

J. Chem. Eng. Data 58, 2904 (2013). 

G. M. Sheldrick, Acta Cryst. A 64, 112 (2008). 

L. J. Farrugia, J. Appl. Cryst. 45, 849 (2012). 

K. Koepernik and H. Eschrig, Phys. Rev. B 59, 1743 
(1999). 

J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992). 
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

V. 1. Anisimov, J. Zaanen, and O. K. Andersen, Phys. 
Rev. B 44, 943 (1991). 

V. 1. Anisimov, F. Aryasetiawan, and A. 1. Lichtenstein, 
J. Phys.: Condens. Matter 9, 767 (1997). 

J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 
118, 8207 (2003); O. A. Vydrov, J. Heyd, A. V. Krukau, 
and G. E. Scuseria, J. Chem. Phys. 125, 074106 (2006). 

'^® G. Kresse and J. Furthmiiller, Comput. Mater. Sci. 6, 15 
(1996); Phys. Rev. B 54, 11169 (1996). 

M. Gentsch and K. Weber, Acta Cryst. C 40, 1309 (1984). 
S. Lebernegg, A. A. Tsirlin, O. Janson, and H. Rosner, 
Phys. Rev. B 89, 165127 (2014). 

S. Lebernegg, M. Schmitt, A. A. Tsirlin, O. Janson, and 

H. Rosner, Phys. Rev. B 87, 155111 (2013). 

H. J. Xiang, E. J. Kan, S.-H. Wei, M.-H. Whangbo, and 

X. G. Gong, Phys. Rev. B 84, 224429 (2011). 

See Supplementary Material for single-crystal crystallo¬ 
graphic data collected at temperatures between 100-280 K, 
the back ground contributions to Cp{T) at H=0T, pow¬ 
der XRD pattern of Langite and exchange constants Jij 
from HSE06 hybrid DFT-calculations, description of the 
tetramer model provided as an alternative to the model 
of two magnetic sublattices, band structures and density 
of states of different structures and optimizations, struc¬ 
tural details of the layers (angles and distances), hydrogen 
bonds and bonding of H 2 O and the SO 4 groups to the 
Cu-O layers;. 

M. Braden, G. Wilkendorf, J. Lorenzana, M. Ain, G. J. 
McIntyre, M. Behruzi, G. Heger, G. Dhalenne, and 
A. Revcolevschi, Phys. Rev. B 54, 1105 (1996). 

A. A. Tsirlin, O. Janson, S. Lebernegg, and H. Rosner, 
Phys. Rev. B 87, 064404 (2013). 

A. Yogi, N. Ahmed, R. Nath, A. A. Tsirlin, S. Kundu, 
A. V. Mahajan, J. Sichelschmidt, B. Roy, and Y. Fu- 
rukawa, Phys. Rev. B 91, 024413 (2015). 


13 


N. Ahmed, A. A. Tsirlin, and R. Nath, Phys. Rev. B 91, 
214413 (2015). 

D. C. Johnston, R. K. Kremer, M. Troyer, X. Wang, 
A. Klumper, S. L. Bud’ko, A. F. Panchula, and P. C. 
Canheld, Phys. Rev. B 61, 9558 (2000). 

For example:Y. C. Arango, E. Vavilova, M. Abdel-Hafiez, 

O. Janson, A. A. Tsirlin, H. Rosner, S.-L. Drechsler, 
M. Weil, G. Nenert, R. Klingeler, O. Volkova, A. Vasiliev, 
V. Kataev, and B. Buchner, Phys. Rev. B 84, 134430 
( 2011 ). 

S. Lebernegg, A. A. Tsirlin, O. Janson, R. Nath, 
J. Sichelschmidt, Y. Skourski, G. Amthauer, and H. Ros¬ 
ner, Phys. Rev. B 84, 174436 (2011). 

F. Illas, M. L, G. de Graaf, and V. Barone, Theor. Ghem. 
Acc. 104, 265 (2000). 

P. J. Hay, J. G. Thibeault, and R. Hoffmann, J. Am. 
Ghem. Soc. 97, 4884 (1975). 

L. Noodleman, J. Ghem. Phys. 74, 5737 (1981). 

A. A. Tsirlin, O. Janson, and H. Rosner, Phys. Rev. B 
84, 144429 (2011). 

C. Yasuda, S. Todo, K. Hukushima, F. Alet, M. Keller, 

M. Troyer, and H. Takayama, Phys. Rev. Lett. 94, 217201 
(2005). 

H.-J. Koo, R. K. Kremer, and M.-H. Whangbo, J. Phys. 
Soc. Jpn. 81, 063704 (2012). 

S. Hara, H. Kondo, and H. Sato, J. Phys. Soc. Jpn. 80, 
043701 (2011). 


14 

Supporting Material 


TABLE III. Obtained lattice parameters of Langite, details of the refinement procedure and settings for the single-crystal 
XRD at various temperatures. 


T(K) 

280 

250 

220 

140 

100 

Crystallographic data 

S.G. 

Pc 

Pc 

Pc 

Pc 

Pc 

a (A) 

7.1412(8) 

7.1370(8) 

7.1349(8) 

7.1255(8) 

7.1231(8) 

6(A) 

6.0468(7) 

6.0427(7) 

6.0400(7) 

6.0312(7) 

6.0305(7) 

C (A) 

11.2328(12) 

11.2238(12) 

11.2204(12) 

11.1996(12) 

11.1935(12) 

P (deg) 

90.0716 

90.0743(14) 

90.0892(14) 

90.1279(14) 

90.1479814) 

Volume (A^) 

485.048 

484.05(9) 

483.54(9) 

481.31(9) 

480.83(9) 

Z 

2 

2 

2 

2 

2 

Dx (mg/m®) 

3.274 

3.281 

3.285 

3.300 

3.303 

Absorption (cm“^) 

8.924 

8.942 

8.951 

8.993 

9.002 

Radiation 

Mo(Kc.) 

Mo(Kc.) 

Mo(Kc) 

Mo(Kc) 

Mo(Kc) 

Data collection 

Theta range (deg) 

2.85-29.00 

2.85-29.01 

2.86-28.87 

2.86-28.97 

2.86-28.97 

Reflection collected 

5811 

5773 

5764 

5726 

5657 

Independent reflections 

2366 

2358 

2359 

2355 

2334 

Completeness of data 

99.8 

99.8 

99.8 

99.8 

99.8 

R^ut (%) 

1.94 

1.95 

2.62 

4.1 

5.2 

Ri (7 > 2a(/)) 

3.28 

3.46 

3.85 

5.2 

4.64 

wR2 

6.33 

6.84 

7.39 

8.01 

6.77 

Largest diff. peak/hole (exA®) 

0.53/-0.75 

0.99/-0.86 

1.12/-0.63 

0.97/-0.90 

1.16/-0.88 


TABLE IV. Refined atomic positions (in fractions of lattice parameters) and isotropic atomic displacement parameters Uiso 
(xlO“^ A^) of Langite collected at 140 K. All atoms are in the general position 2a of the space group Plcl. OW denotes O of 
the H 2 O molecules. OS atoms belong to the SO 4 tetrahedra. 


Atom 

x/a 

y/b 

z/c 

Uiso 

Cul 

0.99970(15) 

0.9975(2) 

0.49967(11) 

0.61(3) 

Cu2 

0.99299(15) 

0.4927(2) 

0.50197(10) 

0.55(3) 

Cu3 

1.00363(16) 

0.7559(2) 

0.75290(12) 

0.53(2) 

Cu4 

1.00899(16) 

0.25469(18) 

0.75168(12) 

0.47(2) 

S 

0.5775(5) 

0.1846(4) 

0.4204(3) 

0.82(4) 

01 

0.8870(10) 

0.9993(11) 

0.6637(7) 

0.68(14) 

02 

0.8874(10) 

0.5052(11) 

0.6661(7) 

0.50(14) 

03 

1.1155(10) 

0.5094(11) 

0.8439(7) 

0.59(15) 

04 

1.1402(9) 

0.2462(12) 

0.5597(7) 

0.48(14) 

05 

0.8583(9) 

0.7436(10) 

0.4428(6) 

0.63(14) 

06 

1.1253(9) 

1.0059(11) 

0.8407(7) 

0.40(15) 

OSl 

0.7837(9) 

0.2316(10) 

0.4119(6) 

0.61(13) 

OS2 

0.5374(9) 

0.0648(11) 

0.5333(7) 

1.38(15) 

OS3 

1.4760(10) 

0.4033(11) 

0.4229(7) 

1.27(15) 

OS4 

1.5146(10) 

0.9525(11) 

0.8181(7) 

1.55(15) 

OWl 

1.2595(10) 

0.7411(11) 

0.6016(8) 

0.94(15) 

OW2 

0.5157(11) 

0.4272(13) 

0.6982(8) 

1.60(17) 
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TABLE V. Refined atomic positions (in fractions of lattice parameters) and isotropic atomic displacement parameters Uiso 
(xlO“^ A^) of Langite collected at 220 K. All atoms are in the general position 2a of the space group Plcl. OW denotes O 
and H atoms of the H 2 O molecules. OS atoms belong to the SO 4 tetrahedra. 


Atom 

xja 

y/b 

z/c 

Uiso 

Cul 

0.99948(9) 

0.99780(11) 

0.49992(6) 

0.926(18) 

Cu2 

0.99321(9) 

0.49226(12) 

0.50203(6) 

0.897(19) 

Cu3 

1.00292(10) 

0.75525(12) 

0.75311(7) 

0.892(15) 

Cu4 

1.00889(10) 

0.25436(11) 

0.75165(7) 

0.830(17) 

S 

0.5780(3) 

0.1835(2) 

0.4194(2) 

1.15(3) 

01 

0.8890(7) 

0.9981(6) 

0.6620(5) 

0.97(9) 

02 

0.8896(7) 

0.5053(6) 

0.6646(4) 

0.88(9) 

03 

1.1182(7) 

0.5106(6) 

0.8425(5) 

0.92(10) 

04 

1.1406(6) 

0.2465(7) 

0.5599(4) 

0.85(10) 

05 

0.8597(7) 

0.7453(6) 

0.4432(4) 

0.98(10) 

06 

1.1249(6) 

1.0072(6) 

0.8388(4) 

0.69(10) 

OSl 

0.7825(7) 

0.2318(6) 

0.4116(4) 

1.13(9) 

OS2 

0.5400(6) 

0.0616(8) 

0.5320(5) 

2.49(11) 

OS3 

1.4768(6) 

0.3973(7) 

0.4229(4) 

1.96(10) 

OS4 

1.5172(7) 

0.9495(8) 

0.8164(5) 

2.44(11) 

OWl 

1.2635(7) 

0.7403(7) 

0.6001(5) 

1.82(10) 

OW2 

0.5190(8) 

0.4213(8) 

0.6976(5) 

2.61(11) 


TABLE VI. Refined atomic positions (in fractions of lattice parameters) and isotropic atomic displacement parameters Uiso 
(xl0“^ A^) of Langite collected at 250 K. All atoms are in the general position 2a of the space group Plcl. OW denotes O 
and H atoms of the H 2 O molecules. OS atoms belong to the SO 4 tetrahedra. 


Atom 

xja 

Vib 

zjc 

Uiso 

Cul 

0.99955(8) 

0.99755(10) 

0.49986(6) 

0.944(16) 

Cu2 

0.99336(8) 

0.49193(10) 

0.50201(5) 

0.878(17) 

Cu3 

1.00281(9) 

0.75523(11) 

0.75317(6) 

0.897(14) 

Cu4 

1.00877(9) 

0.25437(9) 

0.75180(6) 

0.805(15) 

S 

0.5780(3) 

0.1835(2) 

0.4193(2) 

1.24(2) 

01 

0.8878(6) 

0.9990(5) 

0.6620(4) 

0.78(8) 

02 

0 .8886(6) 

0.5053(5) 

0.6637(4) 

0.85(8) 

03 

1.1165(6) 

0.5096(6) 

0.8416(4) 

0.96(9) 

04 

1.1403(6) 

0.2455(6) 

0.5596(4) 

0.87(8) 

05 

0.8608(6) 

0.7446(5) 

0.4427(4) 

0.83(8) 

06 

1.1250(6) 

1.0060(5) 

0.8381(4) 

0.78(9) 

OSl 

0.7835(6) 

0.2322(5) 

0.4113(4) 

1 .00(8) 

OS2 

0.5399(6) 

0.0620(7) 

0.5316(4) 

2.51(9) 

OS3 

1.4780(6) 

0.3970(6) 

0.4231(4) 

2.08(9) 

OS4 

1.5165(6) 

0.9503(7) 

0.8158(4) 

2.49(10) 

OWl 

1.2626(6) 

0.7402(6) 

0.6006(4) 

1.72(9) 

OW2 

0.5193(7) 

0.4239(8) 

0.6992(5) 

2.65(10) 
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TABLE VII. Refined atomic positions (in fractions of lattice parameters) and isotropic atomic displacement parameters Uiso 
(xlO“^ A^) of Langite collected at 280 K. All atoms are in the general position 2a of the space group Plcl. OW denotes O 
and H atoms of the H 2 O molecules. OS atoms belong to the SO 4 tetrahedra. 


Atom 

xja 

Vlb 

zjc 

iso 

Cul 

0.9997(1) 

0.9976(1) 

0.4999(0) 

0.983(14) 

Cu2 

0.9932(1) 

0.4919(1) 

0.5019(0) 

0.968(15) 

Cu3 

1.0028(1) 

0.7552(1) 

0.7531(1) 

0.938(12) 

Cu4 

1.0090(1) 

0.2541(1) 

0.7517(1) 

0.874(14) 

S 

0.5778(3) 

0.1832(2) 

0.4193(2) 

1.35(2) 

01 

0.8878(5) 

0.9988(4) 

0.6631(4) 

0.99(7) 

02 

0.8882(5) 

0.5036(5) 

0.6648(4) 

0.92(7) 

03 

1.1168(5) 

0.5096(5) 

0.8427(4) 

0.94(8) 

04 

1.1407(5) 

0.2452(5) 

0.5600(4) 

1 .01(8) 

05 

0.8608(5) 

0.7446(4) 

0.4433(3) 

0.91(7) 

06 

1.1243(5) 

1.0060(5) 

0.8392(4) 

0.99(8) 

OSl 

0.7832(5) 

0.2308(5) 

0.4109(3) 

1.17(7) 

OS2 

0.5399(5) 

0.0618(6) 

0.5312(4) 

2.63(8) 

OS3 

1.4770(5) 

0.3950(6) 

0.4232(4) 

2.36(8) 

OS4 

1.5171(6) 

0.9498(6) 

0.8156(4) 

2.64(9) 

OWl 

1.2620(6) 

0.7399(6) 

0.6010(4) 

2.03(9) 

OW2 

0.5177(6) 

0.4177(7) 

0.7000(5) 

3.12(10) 
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TABLE VIII. Selected bond lengths (A) for the crystal structnres of Langite collected at different temperatures. 

T(K) 


280 

250 

220 

140 

100 

Cul 

04 

1.927(3) 

1.925(4) 

1.923(4) 

1.924(7) 

1.924(7) 


06 

1.931(3) 

1.931(4) 

1.930(4) 

1.940(6) 

1.938(6) 


oi 

2.000(4) 

1.987(5) 

1.982(6) 

2.007(8) 

2.003(7) 


07 

2.013(4) 

2.024(5) 

2.017(5) 

1.995(8) 

2.006(7) 


08 

2.319(3) 

2.319(4) 

2.319(4) 

2.307(6) 

2.302(6) 

Cu2 

06 

1.914(3) 

1.916(4) 

1.918(4) 

1.910(6) 

1.913(6) 


04 

1.939(3) 

1.933(4) 

1.932(4) 

1.930(7) 

1.940(6) 


02 

1.979(4) 

1.964(5) 

1.970(4) 

1.989(7) 

1.977(7) 


03 

1.995(4) 

2.004(5) 

2 .000(6) 

1.976(8) 

1.983(8) 


08 

2.406(3) 

2.397(4) 

2.401(4) 

2.391(6) 

2.385(6) 

Cu3 

Ol 

1.966(3) 

1.973(4) 

1.964(4) 

1.960(7) 

1.970(6) 


03 

1.970(3) 

1.961(4) 

1.966(5) 

1.970(7) 

1.961(7) 


02 

1.992(3) 

1.988(4) 

1.980(4) 

1.978(7) 

1.984(6) 


07 

1.997(3) 

1.991(4) 

1.999(4) 

1.997(7) 

1.990(7) 


08 

2.368(4) 

2.368(5) 

2.375(5) 

2.376(7) 

2.370(6) 

Cu4 

07 

1.974(4) 

1.969(4) 

1.967(4) 

1.982(7) 

1.980(7) 


02 

1.993(4) 

2.003(4) 

1.994(4) 

1.987(7) 

1.984(6) 


03 

2.006(4) 

1.996(4) 

2.011(5) 

2.000(7) 

1.996(7) 


Ol 

2.031(3) 

2.035(4) 

2.035(5) 

2.024(7) 

2.025(6) 


04 

2.351(4) 

2.353(4) 

2.348(5) 

2.347(8) 

2.344(7) 


06 

2.398(4) 

2.389(5) 

2.398(5) 

2.397(7) 

2.372(7) 

SI 

on 

1.470(4) 

1.475(4) 

1.480(5) 

1.506(7) 

1.488(63) 


012 

1.480(5) 

1.482(5) 

1.473(5) 

1.482(8) 

1.491(7) 


OlO 

1.481(5) 

1.484(5) 

1.487(6) 

1.485(8) 

1.493(7) 


08 

1.498(4) 

1.499(5) 

1.491(5) 

1.500(7) 

1.501(6) 


TABLE IX. The total nearest-neighbor Jij (K) calcnlated with the HSE06 hybrid functional as implemented in VASP5.2 in 
comparison with the results from LSDA-l-17 calculations using the f plo9.07-41 code, Ud = 8.5 ± 1 eV and Jd = 1 eV. d(Cu-Cu) 
and Cu-O-Cu denote Cu-Cu distances (A) and Cu-O-Cu angles (deg), respectively. 


d(Cu-Cu) 

Cu-O-Cu 

HSE06 

LSDA-fl7 

Js 

2.983 

101.05 

-116 

-74± 10 

J[ 

3.011 

97.81/98.64 

21 

9± 15 

Ji 

3.020 

99.10/99.49 

58 

38 ±20 

J's 

3.049 

104.706 

-65 

-23 ±3 

Ja 

3.139 

103.55 

-52 

-19 ± 1 

Jb 

3.144 

105.08 

-23 

5±5 

Jc 

3.163 

106.64 

-43 

-6 ±2 

Jd 

3.168 

105.23 

-16 

8±5 

Je 

3.186 

106.65 

-33 

0±4 

Jf 

3.190 

106.61 

-37 

-12 ±2 

Jq 

3.219 

106.10 

-37 

-13 ± 1 

Jh 

3.229 

109.23 

5 

31 ±7 
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Temperature (K) 

FIG. 7. (Color online) Specific heat data C'p(T)of Langite collected in zero magnetic field. The blue line shows the lattice 

Tl —7 

background CiatiT) which we obtained by fitting a polynomial of the form Ciat{T) — X] to the experimental data in the 

n=3 

temperature regime 20-40K. Cp{T) — Ciat{T) yields the magnetic contribution Cmag to the specific heat. 
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Experimental pattern 
[96-900-9716] Langite 
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FIG. 9. (Color online) The LDA band structure (orange) computed for the crystal structure provided by Gentsch and Weber 
in 1984 (GW). Their data was collected at room temperature (RT). The bands computed for GGA-optimized H-positions (opt 
H), with all other crystallographic data being fixed to those provided by GW, are shown in green. The optimization stabilizes 
bands up to about 50 meV and significantly changes the band dispersion which is particularly pronounced in the energy region 
between -0.5 and -0.2 eV. 
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FIG. 10. (Color online) The LDA band structure (green) computed for the crystal structure provided by Gentsch and Weber 
in 1984 (GW) with optimized H-positions (opt H). The red bands are obtained from our crystallographic data with optimized 
H-positions (opt H). In both cases one H is located at the SO 4 group. The differences probably arise from the different 
temperatures. While GW collected their data at room temperature (RT) ours was measured at 100 K. 
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FIG. 11. (Color online) The LDA band structure (red) computed for our crystal structure with optimized H-positions (opt 
H). The blue bands correspond to the structure where H and the oxygen atoms of the water molecules (OW) are optimized. 
In the hrst case one H is located at the SO 4 group. While in the latter case (opt H+OW) H from the SO 4 group has moved 
to the Cu-O layers. Particularly in the energy regime between -0.6 to -0.2 eV band shifts up to 100 meV are visible as well as 
considerable changes in the band dispersion. The structure with no H on SO 4 is considerably more stable. 














































23 



Temperature (K) 


FIG. 12. (Color online) A comparison of the density of states around the Fermi level calculated within the LDA for different 
crystal structures. ”GW”: the structure refined by Gentsch and Weber in 1984 at room temperature. ”GW opt H”: The GW 
crystal structure with GGA-optimized H-positions. ’’opt H” denotes our crystallographic data collected at 100K with optimize 
H-positions and ’’opt H+OW” that with optimized H and OW positions, where OW is the oxygen of the water molecules. 



b 


FIG. 13. (Color online) The bonding of the H 2 O molecules and the SO4 groups. Thick and thin dashed blue lines denote 
H-bonds shorter and longer than 0.80 A. The dashed brown lines indicate the weak bonds between O atoms belonging to H 2 O 
or SO4 and Cu atoms. 
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FIG. 14. (Color online) The right figure shows the structural layer in langite with the four different types of Cu 04 plaquettes. 
Black lines indicate Cu-O bonds. The black numbers with three digits denote the corresponding Cu-O bond lengths (in A). 
Blue numbers with two digits give the Cu-O-Cu bond angles in deg. To the left, fragments of type-A and type-B chains are 
shown, where the numbers give the O-Cu-0 angles (in deg) on the different plaquettes. 
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FIG. 15. (Color online) The left panel shows the intralayer exchange pathways. Transparent green areas represent the positions 
of the 5" = 1 tetramers. These tetramers and the belonging nearest-neighbor are shown in the central graphic where solid 
and dashed lines indicate AFM and FM interactions, respectively. J/ is depicted in light green color since this coupling is 
frustrated. The weak couplings Jj, and Ja are not displayed. Red arrows indicate the alignment of the Sz = 1/2 spins for the 
isotropic model. The right panel shows the two chain types and the respective intrachain exchange pathways. 

Description of the tetramer model 

Considering first only the nearest-neighbor couplings (see Fig. 15), an almost unfrustrated microscopic magnetic 
model for langite is obtained. Only the weak couplings Jb and J/ cannot be satisfied. Jg and J' are both FM 
with J/ being much smaller than Jg rendering the B-type chains alternating ferromagnetic in this simple model. 

In the edge-sharing type A chains Ji is AFM and about four times larger than the weak J[. Thus, the A-type 
chains can be reduced to AFM Cu3-Cu4 dimers. Ju, reigning between Cu2 and Cu3, couples these dimers 
to the ferromagnetic B-type chains resulting in S' = 1 tetramers (central panel in Fig.S3). Cul and Cu3 are 
coupled ferromagnetically by Ja resulting in an antiparallel alignment of the tetramers along c. J', in turn, 
is responsible for a parallel alignment along the chain direction b. Now, next-nearest neighbor exchanges are 
taken into account: J/ and J// can be satisfied within the tetramer model but J 2 and J 2 are frustrated. The 
model of coupled S = 1 tetramers allows the following phenomenological interpretation of x{T) and M(H). The 
bend in M(H), starting at about 10 T, arises from a parallel alignment of the S = 1 tetramers in the magnetic 
field. In this case, 3/4 of the spins are oriented parallel to the magnetic field resulting in a magnetization 
of O.b/rs/Cu in agreement with the measurements. The effective AFM coupling between the tetramers Ju 
can roughly be estimated by an AFM S' = 1 square-lattice spin model yielding SJukB/~ lOT and 
accordingly Ju « 1.8 K. This effective coupling is considerably smaller than our calculated values, e.g. for J/ 
and J//, however, the seeming discrepancy can be explained by frustration resulting in the very weak effective 
intertetramer coupling. A similarly small value for Ju can be obtained from x(T). The maximum of the 
magnetic susceptibility for an AFM S = 1 square-lattice is at about 2.1%T/J^. According to the experimental 
x(^)max = 7.5 K, we, thus, obtain Ju = 3.4 K. With respect to our crude estimate based on an AFM S = 1 
square-lattice the agreement between the Ju extracted from x{T) and M(H) is reasonable and, therefore, justifies 
the S' = 1 tetramer model and also clearly shows the crucial role of frustration in langite. 

We employ again the tetramer model and assume an intratetramer coupling of about 20-40 K. Now, a saturation 
field may roughly be estimated using a two-dimensional system with four effective couplings on each Cu-site 
yielding ^Jks/{gg^B) « 60 — 120T. The bend at about 4T cannot easily be ascribed to certain mechanisms 
within the phenomenological model but may arise from secondary effects. 

The model of coupled tetramers cannot provide a better description of the magnetic properties than the model 
of coupled magnetic sublattices which we present in the main text. It also neglects the strong NNN J 2 and J 2 
couplings and appears to be counterintuitive since chain features are completely lost. Such feature, however, were 
observed in neutron experiments on related compounds like brochantite and antlerite (see main text). Therefore, 
it appears natural to stick to a model based on magnetic chains as long as it is not proved inappropriate. 

















